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Abstract 

In lizards, one of the most important behavioral mechanisms to cope with spatial and temporal variations in thermal 
resources observed is activity time. The longer a lizard can maintain activity, the more time it has to forage and reach larger 
adult body size. We studied the behavioral adjustments to different climatic regimens on daily and seasonal scales in three 
natural populations of the lizard Microlophus atacamensis along a latitudinal temperature and rainfall gradient. We also used 
Niche Mapper to determinate the amount of thermally suitable time for activity for this species. Abundance and daily 
activity patterns varied greatly over the year for the three populations. In summer and spring, the daily activity times were 
greater, and were reduced in fall and winter seasons. In summer, when stressful heat loads should prohibit activity over a 
midday gap, lizards did not show bimodal patterns of activity. Instead, they move to the cooler intertidal habitat. 
Abundance and thermal quality in the southernmost coolest site was lower, and the potential annual activity time decreases 
with latitude. Contrary to expectations, lizards from this locality showed the largest body sizes possibly due to diet and/or 
time to sexual maturation. Our results indicate that the intertidal habitat is a key factor that influences daily and seasonal 
activity of M. atacamensis lizards. While this habitat is not climatically optimal for lizards, it allows them to behaviorally 
extend their activity window and gain access to food in the intertidal areas. 
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Introduction 

One of the principal ciiallenges for animals is to obtain enough 
food to provide nutrients and energy needed for maintenance, 
survival, growth, and reproduction [1]. The energy allocated to 
these processes ultimately depends on tlie daily and annual energy 
budget, which in turn depends on activity time [2]. In lizards, it 
has been demonstrated that the daily and seasonal activity times 
may be influenced by both endogenous and exogenous environ- 
mental factors [3] , [4] . The endogenous rhythms of activity are 
driven both by circadian and circannual clocks that allow changes 
in the duration of daily activity across seasons, allowing the lizard 
to anticipate the upcoming season [5]. However, environmental 
exogenous factors may mask endogenous rhythms by either 
increasing or suppressing activity, the most important being 
ambient temperatures in the field, which influence both daily and 
seasonal activity patterns [3], [6], [7], [8]. Assuming that lizards 
are foraging while active, and that food is not limiting in a given 
environment, changes in daily or seasonal activity times due to 
spatial and temporal variations in thermal resources influence the 
daily prey capture rate [2], [9], [10], thus ultimately affecting 
energy acquisition and the reproductive success of lizards. 



Daily and seasonal activity may also lead to specific geographic 
patterns in growth and body size. Lizards with longer activity 
seasons are expected to grow faster [2]. For example, in sagebrush 
lizards, it was found that longer seasonal activity was associated 
with an increase in individual growth and larger adult body size 
[1 1] (but see [2]). Therefore, variation in thermal environment, 
through its effect on thermoregulatory behavior and activity time, 
may have a profound and direct impact on time budgets of the 
lizards, and cause some of the observed variation in life histories 
among widespread populations [2], [9], [12], [13]. 

The family Tropiduridae is mostly composed by terrestrial feeders, 
but some species that dwell in desert coastal areas (with low levels 
of terrestrial productivity), feed almost exclusively on intertidal 
zones [14], [15]. In Chile, two species from the genus Microlophus 
use this feeding strategy: Microlophus quadrivittatus and M. atacamensis 
[16]. The geographical range of M. quadrivitattus (from Arica 
18°00'S to Antofagasta 20°32'S) falls into an area where the 
extreme desert reaches the coast and there is practically no annual 
precipitation (maximum recorded =2 mm year" ) [17], whereas 
the geographical range of M. atacamensis, the subject of this study 
and southernmost species of the genus, runs 500 km into the so- 
called Desert Coastal area which presents a gradient of 
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precipitation from 2 mm year ' in the North (at Antofagasta) to 
almost 120 mm year"' in the South (La Serena 29°55'S) [18]. For 
this species it is reported that the diet of this species is composed 
both of intertidal (algae and Crustacea) and terrestrial (isopoda, 
coleoptera) prey [18]. The importance of each of these dietary 
items is associated with spatial and temporal variation in the 
productivity levels of both habitats, with an increasing proportion 
of terrestrial items in a north to south gradient [18], [19]. 

The intertidal zone is an unfavorable thermal environment 
because lizards lose heat quickly by conduction, convection and 
evaporation when feeding on a wet rock substrate [20]. For that 
reason, lizards may face a trade-off between feeding and 
thermoregulation. Because M. atacamensis populations are under 
different thermal constraints at latitudinal, seasonal and micro- 
habitat scales, a temporal and geographical variability in their 
thermal biology is expected. Lack of physiological intraspecific 
variation in thermal biology in M. atacamensis (e.g. selected body 
temperatures and heating and cooling rates) suggests that lizards 
may compensate adequately the variations in their thermal 
environment through a behavioral mechanism instead of physi- 
ological adjustment [21]. 

Through field observations and using a process-driven mech- 
anistic model (Niche Mapper) [22] we explored the behavioral 
adjustments of M. atacamensis to diflFerent climatic regimens, and 
their potential for activity at different temporal and spatial scales. 
We selected three populations along the entire geographical range 
of this species. Considering that environmental temperature 
decreases as latitude increases, we specifi[:ally assessed whether 
daily and seasonal activity patterns are modified with increasing 
latitude. Associated with this hypothesis, it is expected that body 
size should decrease in populations that exhibit shorter activity 
times. 

Materials and Methods 

Organism and study sites 

From July 2005 to April 2006 we studied the abundance and 
activity patterns in three sites on the northern Chilean coast: 
Medano (24°37'S; 70°33'W), Zenteno (26°5rS; 70°49'W), and 
Arrayan (29°4rS; 71°19'W). Climatological data for the study 
sites are presented in Table 1. Medano and Zenteno are under a 
coastal desert climatic regime, characterized by less than 2 mm of 
rain per year [17]. The southernmost locality, Arrayan, in 
contrast, falls into the per-arid Mediterranean zone. Rainfall in 
this site averages 120 mm/year. Terrestrial productivity also 
shows a gradient, with a significant increment of plant cover, from 
ca. 0% in Medano to 2'i"/« in Arrayan [18]. Lizards live along the 
intertidal zone in a stretch of approximately 150-200 m domi- 
nated by flat areas of a mixture of rocks (with heights from 0.2 to 
3 m), pebbles, cobblestones, and sand [23]. 

Abundance 

For each site, and every four months, the abundance of active 
M. atacamensis lizards at intertidal and terrestrial areas was assessed 
by visual counts. On each occasion, and for two consecutive days, 
censuses data were collected by walking parallel to the shoreline, 

covering both intertidal and the adjacent terrestrial (supralittoral) 
substrates. Censuses were made at 0900, 1100, 1300, 1500, and 
1700 h, walking on each occasion for one hour and covering 
approximately 1 km of coast. The maximum number of lizards 
sighted during a census period was used as a relative index of 
abundance. This index does not provide a measure of the actual 
density of lizards, but allows the between-habitats and between- 
seasons comparison of abundance at a geographical scale [24]. 



Lizards were differentiated by age class (juveniles or adults), 
recording their position on intertidal or terrestrial habitats. Air and 
soil temperatures (T^ and Tsoii) were measured immediately prior 
and after each census at the beginning and end of the transect. 

Dally activity patterns 

Simultaneously with the abundance analysis, one observer 
counted lizards from 0900 to 1800 h in a 1500 m^ area, recording 
the number of lizards sighted ever^' 30 min. We selected this time 
window because M. atacamensis is a diurnal lizard that restricts their 
activity to daytime hours when Ts„;i is above 22°C [18]. As above, 
lizards were differentiated according to age classes (juveniles or 
adults) and substiate (intertidal or terrestrial habitat). To establish 
an approximation of thermal conditions of the habitat occupied by 
lizards, we placed six Tidbit temperature data loggers with an 
internal temperature sensor (three in the intertidal zone and three 
in terrestrial zones) at randomly selected microsites for each season 
and locality. Temperatures in the data loggers were registered 
every 30 min. Thermal profiles collected from the electronic 
devices do not differ substantially from copper lizards models 
exposed to the same conditions, suggesting that Tidbit could be 
used directiy to monitor operative temperatures (TJ [25]. To test 
if this assumption was true for our model species, in the laboratory 
we calibrated the tidbits measurements with those of dead lizards 
placed in a thermal gradient. The results showed that data 
recorded from electronic devices are highly correlated with 
temperatures from dead M. atacamensis lizards exposed to identical 
thermal conditions (r= 0.99).We also calculated the d,, index, which 
represents the thermal quality of the habitat, from the lizard's 
perspective [26]. This index is calculated as the mean of the 
absolute value of the deviations of Tg from selected temperature 
(Tjpi) of the corresponding site and season. It was found that M. 
atacamensis lizards have a T^,.! of 34.3°C, independent of localities 
or seasons [21]. For this reason, we use this value for thermal 
quality measurements. 

Activity times 

To estimate the maximum duration of daily and annual activity 
of lizards from the three sites, we used the Niche Mapper model 
[6], [22], to determinate the biophysical eflFe[:ts on the amount of 
thermally suitable time for activity, under realistic climatic 
conditions [27]. The microclimate model in Niche Mapper allows 
the translation from macro to microclimatic environmental 
variables relevant to the organisms. More details of these models 
can be found in [6], [22], and at [http://www.zoology.wisc.edu/ 
faculty/Por/Por.html] . 

The microclimate model contains a subroutine for l onijjuting 
solar radiation given a specific date, time, latitude, longitude, 
slope, and aspect [28]. Time dependent input data consist of 
monthly 2 m shade air temperatures, wind speed, and relative 
humidities. These data for the three sites were obtained from the 
Servicio Hidrologico y Oceanografico of the Armada de Chile 
(Table 1). The microclimate model calculates the hourly intensity 
of solar radiation and the profiles of and Tsoii! wind and 
humidity above the surface. These data were used for calculating 
the hottest (maximum sunlight available) and coolest (maximum 
shade available) microclimates available to the lizards in the rocks 
for each hour of the day and for each location. Hourly output from 
the microclimate model for the hottest and coolest above and 
below ground conditions were used to define the fuU range of 
environmental conditions available locally. Time dependent input 
data used for the microclimate simulations can be found in 
Tables 1, 2 and 3. The microclimate output conditions were then 
used as input for the ectotherm model, in addition to various 
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Table 1. Climatological data for the study sites in the four seasons of the year. 





Season 




Locality 










Medano 


Zenteno 


Arrayan 


Summer 


Tmax 


22.9±0.8 


23.2 ±1.7 


20.0 ±1.8 






18.6±1.0 


16.9±1.2 


15.8±0.9 






20.7±1.7 


19.6±2.5 


17.6±1.8 


Fall 


iir 


19.7±2.3 


20.0±2.3 


16.9 ±2.4 




Tmin 


15.6±2.3 


14.0±2.4 


12.7±2.4 




Tmean 


17.6±2.7 


16.8±2.9 


14.5 ±2.6 


Winter 


tII" 


16.5±1.1 


17.9±1.5 


15.4±1.5 




Tmin 


13.2±1.1 


11.9±2.1 


11.1 ±1.4 




Tmean 


14.9±1.5 


14.6±2.4 


13.0±1.9 


Spring 


Tmax 


18.9±1.6 


19.9±1.9 


17.4±1.6 




Tmin 


16.0±1.3 


14.1 ±2.4 


13.3±1.3 




Tmean 


17.5±1.8 


16.7 ±2.3 


15.0±1.9 


Relative humidity {%) 




85.6 




84.5 


Precipitation (mm) 




2.2 


2.2 


127.4 


Climate 




Desert 


Desert 


Mediterranean 



Note: Values are expressed as mean ± SD. Data were obtained from the Servicio Hidrologico y Oceanografico of the Armada de Chile. Only precipitation was obtained 
from [17]. 

doi:l 0.1 371 /journal.pone.0097735.t001 



properties of the animal relating to morphological, physical, 
physiological, food, and behavioral properties [27]. 

The ectotherm program uses the results from the microclimate 
model to determine which hours of the day an ectotherm can stay 
within its preferred body temperature range. We used this model 
to examine the frequency distribution of expected daytime body 
temperature (T|,) of M. atacamensis lizards throughout the year in 
the three sites. The model assumes that the animal is always active 
(and exposed to local wind and humidity) when thermal conditions 



permit. If conditions are too hot or cold for activity, the animal 
retreats (close to the sea or underground) to a temperature as close 
to its optimal T|, as possible. Analyses were done separately for wet 
and dry conditions, simulating both intertidal and terrestrial 
substrates, respectively. Calculations for the wet substrates 
assumed a 100% wet surface. Wet and dry reflectivities of local 
substrates were measured using an ASD portable spectroreflect- 
ometer for the energy balance calculations for all habitats. The 
spectral range of these measurements was 350-2500 nm, approx- 



Table 2. Parameters used in Niche IVlapper for Microlophus atacamensis lizards. 



Variable 


Value 


Source 


Mean male/female weight, Medano (g) 


65.2/35.7 


This study 


Mean male/female weight, Zenteno (g) 


81.2/35.8 


This study 


Mean male/female weight, Arrayan (g) 


84.9/45.8 


This study 


Oxygen extraction efficiency {%) 


6-12 


Assumed 


Activity above basal multiplier 


3 


Assumed 


Preferred body temperature, C 


33.7-34.7 


[21] 


Voluntary thermal minimum, C 


28 


[21] 


Voluntary thermal maximum, C 


37 


[21] 


Fecal water (%) 


10 


Assumed 


Dry matter of food (%} 


31 


Assumed 


Dry mass of food carbohydrate (%} 


12 


Assumed 


Dry mass of food lipids (%) 


18 


Assumed 


Dry mass of food protein {%) 


67 


Assumed 


Activity period 


Diurnal 


This study 


Percent skin wet 


0 


Assumed 


Maximum/minimum absorptivity (%) 


92/85 


ASD measurements. This study 



doi:1 0.1 371 /journal.pone.0097735.t002 
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imately 97% of the sunlight reaching the surface of the earth. Rock 
reflectivities were 21/32% (dark/light) for Medano, 21% for 
Zenteno, and 16% for Arrayan. Input data for the ectotherm 
simulations can be found in Table 2. 

Body size 

Live lizards from the three populations were captured during 
the study period (N = 55 for Medano, N = 50 for Zenteno, and 
N 27 for Arrayan). Two measures of body size were registered at 
each capture: (1) snout-vent length (SVL), which was measured to 
the nearest millimeter using a caliper, and (2) mass, which was 
measured to the nearest 0.1 g using a Pesola scale. The analyses 
were restricted to adults and, because this species exhibits a clear 
sexual dimorphism in body size [16], [23], sexes were analyzed 
separately. As soon as the measurements were taken the 
individuals were liberated. 

Ethics statement 

The protocol of capture was done with research permits from 
the Chilean Agriculture and Livestock Bureau (SAG) (Permit 
Number 98), and under the approval of the Pontificia Universidad 
Catolica de Santiago Ethics Committee. No specific permissions 
were required in the study sites because they are within an open 
access zone of the Atacama Desert coast. Also, M. atacamensis is 
classified as Vulnerable under the Chilean Hunting Law (DS 5). 

Statistical analyses 

The abundance of At. atacamensis was analyzed with a two-way 
ANOVA [29], using sites and seasons as factors, with 3 and 4 
levels, respectively. KendaU-rank correlations were used to 
determine whether the mean number of lizards was correlated 
with Tair or Tjoji. Means of T^ir and T^so'i were used to compute 
this statistic. Thermal conditions during the daily activity patterns 
(using the dg index) were compared among sites and seasons using 
a two-way ANOVA. Body size data (SVL and mass) was analyzed 
with a two-way ANOVA, using sites and sex as factors. Normality 
of the data was checked by Shapiro-Wilk's test of normality and, 
when necessary, the data were log-transformed. Post-hoc com- 
parisons were performed using a Tukey's honest significance 
difference test. AU comparisons were performed using Statistica 
8.0 for Windows [30]. 

Results 

Abundance 

The abundance of M. atacamensis was significantly different 
among sites (^"2,12 = 14.40, P<0.001), seasons (i?3,i2 = 9-73, 
P= 0.002), and for the interactions between these factors 
(^6,12 ~ 4.44, P= 0.014). Highest mean values occurred in Medano 
and Zenteno (at the north and mid geographical range, 
respectively) compared to the southernmost site (Arrayan) (Fig. 1). 
In terms of seasonality, fewer individuals were found in winter, 
compared to the summer. T^q^, and T,. were higher in 

summer, related to the other seasons (Table 3).Tjjr and Ts„ij were 
significantiy correlated with the number of lizards observed during 
the censuses, for the three sites studied (t = 0.49 and t = 0.35 in 
Medano; t = 0.35 and t = 0.40 in Zenteno; t = 0.18 and t = 0.14 
in Arrayan, for Tj,;^ and T^,,;!, respectively) (all Kendall rank 
correlation coefficients significant at /'<0.001). 

Daily and annual activity patterns 

Daily activity patterns varied gready among seasons for the 
three populations. During summer, lizards emerged from their 
burrows as early as 0900 h, when T^ was above 20°C. Lizards 
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Fig. 1. Seasonal and geographical comparisons of mean lizard abundance. Vertical lines represent standard errors. 
dor:1 0.1 371/journal.pone.0097735.g001 



remain active all day, and did not show a clear abundance pattern 
(Fig. 2). When we finished our observation period (1800 h), we 
could stiU observe some individuals above the rocks. Both juveniles 
and adults emerged at 0900 h, and remained active until 1800 h. 

The activit)' time was shorter in fall and winter, especially for 
Medano and Arrayan sites. Lizards' appearance coincided with an 
increase in T^. (aljove 25 C), showing a unimodal activity pattern. 
In general, juveniles emerged earlier than adults in the three sites 
and disappear before them in late afternoon. In spring, daily 
activity pattern was similar to fall and winter, showing a unimodal 
activity pattern, with maximal abundances near midday for both 
juveniles and adults. The activity period in spring months was 
similar to summer, but hzards emerged around 1000 h for the 
three sites. Again, juveniles emerged earlier than adults, but 
disappear before them in the afternoon. 



Operative temperature was strongly correlated with the number 
of observed ac:tive hzards, for the three sites studied (t = 0.62 in 
Medano; x — 0.34 in Zenteno; T = 0.57 in Arrayan) (aJl Kendall 
rank correlation coefficients significant at /'<0.001, seasons 
pooled). The deviation of T,. from T,.,.! (d,) varied significantly 
among sites (i^p_2i6 = 9. 19, P<0.001), seasons (F^.zie = 35.52, P< 
0.001), and for the interaction of these two factors (/^6,2i6 = 8.80, 
P<0.001). Thus in Medano, average dg was smallest (i.e. a highest 
thermal quality) than in the other sites. Seasonally, dg index was 
lowest in summer and highest in winter. No differences were found 
between spring and fall (Fig. 3). T^ varied considerably with time 
of day for each season and site. In summer, and for all sites, 
thermal suitability for lizards was high during the day. In contrast, 
in the other seasons thermal suitability was lower, especially in 
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Fig. 2. Number of juveniles (white bars) and adults (blacl< bars) M. atacamensis lizards along the day in summer, fall, winter, and 
spring seasons for: Medano (a-d), Zenteno (e-h), and Arrayan (i-l) populations. Circles represent Te (measured with Tidbits data loggers). 
doi:1 0.1 371/journal.pone.0097735.g002 
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Fig. 3. Mean and standard error of the index. 
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early morning and late aftemoon, and higher at midday hours, 
when lizards were active. 

Microlophus atacamensis also shows geographical and seasonal 
differences in the use of intertidal and terrestrial habitats. In 
Medano, lizards show a clear pattern of use of both habitats. 
Individuals were found in the intertidal zone when Te surpasses 
25°C. In winter, when almost never exceeds this temperature, 
Uzards remain in the terrestrial arc'as (Fig. 4). In summer and 
spring, Kzards from Zenteno and Arrayan were found both in the 
intertidal and terrestrial habitats during the day. In fall, only a few 
individuals were found in the intertidal zone during a short period 
of time. In winter, no individuals were found in the intertidal zone 
for any site. 



Biophysical simulations 

Calculations of core temperature throughout a typical day in 
different months and hours of day were made for two different 
body weights (100 g in males and 36 g in females) in each site, 
which represent the average weights for both sexes' measured 
levels of activity. Considering that the voluntary thermal minimum 
and maximum for lizards are 28 and 37°C, respectively (Table 2), 
constraints on daily and seasonal potential activity times under the 
dry (terrestrial) and wet (intertidal) conditions were estimated 
(Fig. 5). As was observed in the field, our simulations largely reflect 
seasonal variations in activity times in all sites, with winter activity 
significantiy curtailed by temperature. The thermal constraints on 
lizard activity precluded the occurrence of M. atacamensis in the 




MO lOJO VIM 13:30 ISM 1«30 UM MO lOdO U:30 U:00 ItM IMO MO lOrJO 12.-00 13-JO ISM IC-JO IMO 



Haurofdif 

Fig. 4. Diel variation of the percentage of active lizards in terrestrial and intertidal habitats in summer, fall, winter, and spring 
seasons for: Medano (a-d), Zenteno (e-h), and Arrayan (i-l) populations. White bars represent the proportion of individuals In the terrestrial 
or supralittoral zone; the black ones represent the proportion in the intertidal zone. Maximal abundance In the day was considered as 100%. Circles 
represent Te (measured with Tidbits data loggers) in the intertidal and terrestrial habitats. 
doi:1 0.1 371/journal.pone.0097735.g004 
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Fig. 5. Seasonal and daily activity constraints for Microlophus atacamensis'in IVIedano, Zenteno, and Arrayan, on dry (terrestrial) and 
wet (intertidal) conditions. Gray contours indicate suitable temperatures (28-36°C) for activity. 
doi:1 0.1 371 /journal.pone.0097735.g005 



intertidal liabitat in late fall and early winter months, especially for 
Arrayan site. 

Based on the biophysical model, the activity hours for lizards in 
the three localities are highest in spring and summer months, when 
the ability to exploit both habitats is enhanced (Fig. 6). However, 
the activity decreases in fall and winter months, especially in the 
intertidal habitat. Potential annual activity time decreases with 
latitude, ranging from 3099 h yr 'in Medano to 2981 h yr 'in 
Arrayan in dry conditions, and also in wet conditions from 2918 to 
2860 h yr ' in Medano and Arrayan, respectively. Limited 
activity is manifested in reduced opportunities for feeding, which 
can be assessed by measuring tlu' required foraging rates in each 
locality and habitat. In northern sites, and in dry conditions, the 
foraging rate increases because of increased opportunities for 
activity (Fig. 7). In wet and southern habitats, lizards are 
compromised in terms of total annual activity hours, but their 
annual food requirements are reduced. Through the year, foraging 
rates are generally higher (~0. 1 g h~ ') during spring and summer 
months and lower (<0.07 g h ') during fall and winter seasons. 

Body size 

Body size (SVL and mass) differed significantly among sites and 
sexes (Table 4). For SVL comparisons, females were significantiy 
smaller than males in all the sites studied. However, for each sex, 
Uzards from Medano, Zenteno, and Arrayan were similar in size. 
With regard to body mass, males were heavier than females. In 
females, body mass of Arrayan lizards was greater than both 
Medano and Zenteno lizards. Medano and Zenteno lizards were 
similar in weight. In the case of males, body mass from Medano 
lizards was smaller than lizards from Zenteno and Arrayan, but 
similar between Zenteno and Arrayan. 

Discussion 

Our results indicate that the latitude and climate are important 

factors in modulating the activity patterns of M. atacamensis lizards. 
We predicted that hzards from northern sites should have longer 
daily and seasonal activity times because they experience a warmer 



environment. Accordingly, we found a small (4 and a higher 
potential annual activity time in Medano, implying that this 
northern locality offers better quality thermal conditions for the 
activity of lizards [13]. A decrease in potential activity time in the 
southernmost site (Arrayan) is explained because the thermal 
quality of the environment in this locality is lower. It is likely that 
lizards from this site compensate for the high time cost of 
thermoregulation by allocating more time to basking than to other 
activities [31], which in turn may imply a trade-off between 
feeding and thermoregulation. An increase in basking time 
activities can affect foraging time and consequently the energy 
budget [32], and ultimately the reproductive success of lizards in 
high-latitude populations. The abundance of lizards in Arrayan 
was significandy lower than in Medano and Zenteno, suggesting 
that in this southern locality the animals were under a higher cost 
of thermoregulation. Lizards from Arrayan should allocate more 
time to thermoregulation with a consequendy smaller proportion 
of the time budget allocated to foraging, thus decreasing the 
energy available for growth and reproduction. It is important to 
note that the smaller population size in Arrayan compared to 
Medano and Zentano might be also due to higher predation (e.g. 
more birds and snakes in a more mesic environment) on adult and 
juvenile lizards [18]. 

We found that these lizards also shift their daily activity times 
during the year. In summer and spring the daily activity time was 
longer, and reduced in the fall and winter seasons as expected for 
diurnal lizards [6], [33], [34]. Our biophysical model also revealed 
similar results, with a longer period of activity in summer and 
spring and lower in fall and winter, providing good correspon- 
dence with empirical data. Similar to the differences among 
localities, longer daily activity times found in summer and spring 
could be explained because the d,, index was lower in these seasons, 
indicating a high thermal quality for the activity of lizards. 

During the fall, winter and spring seasons, lizards showed a 
unimodal activity pattern, being active around midday when solar 
radiation and T^ were high. This pattern is typical of several 
species of lizards [35] , [36] . In particular, M. pemvianus (a species 




Fig. 6. Monthly predicted hours of activity for lizards in IVIedano, Zenteno, and Arrayan, on dry (terrestrial, filled symbols) and wet 
(intertidal, empty symbols) conditions. 

doi:1 0.1 371/journal.pone.0097735.g006 
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Fig. 7. Foraging rates for lizards from Medano, Zenteno, and Arrayan on dry (terrestrial, filled symbols) and wet (intertidal, empty 
symbols) conditions. 

doi:1 0.1 371 /journal.pone.0097735.g007 



that belongs to the M. pemvianus group, as does M. atacamensis) also 
shows a unimodal pattern of activity during this period of time 
[20]. Nevertheless, a shift in activity from unimodal to bimodal has 
been found in several species of diurnal desert lizards in summer 
[6], [37], [38], generally observed at sites where air and soil 
temperatures reached risky highest levels. Around midday, lizards 
retreat into the shade or burrow to avoid overheating, but become 
active again during the afternoon and thus exhibit a markedly 
bimodal activity pattern [37]. Traditionally, this daily pattern of 
activity has been explained as the behavioral response of these 
animals to the very high levels of air and soil temperatures at that 
time of day [39]. However, laboratory experiments in some species 
of lizards demonstrate that summer bimodality is not merely a 
proximate response of lizards to high levels of soil temperature, but 
is a clock driven phenomenon mediated by circadian rhythmic 
changes of plasma melatonin levels [40], [41], [42], which has 
most probably evolved as a behavioral adaptation to high levels of 
air temperature occurring predictably around midday [5], [43]. 

Regardless of the mechanism associated with a unimodal/ 
bimodal activity pattern, during summer M. atacamensis faces a 
stressful heat load in midday hours, which should prohibit activity 
during these hours. However, they did not show a bimodal pattern 
of activity, but interestingly they used the behavioral strategy of 
moving to the intertidal habitat, where temperatures are lower and 



wind speed on the shore is greater. This strategy of lizards of 
moving to the intertidal habitat as the temperature rises allows 
them to engage in continuous daily activity associated with the use 
of different substrates, extending their activity period close to a 
source of food, maintaining a similar temperature during the day 
[44]. Also, this behavior benefits lizards because food and energy 
requirements in wet habitats are reduced. However, it also 
represents costs for thermoregulation, because the low substrate 
and air temperatures in the intertidal habitat increase the heat lost 
by conduction, substrate radiation and local air temperature due 
to local substrate evaporation [20]. 

Body temperature showed significant variation throughout the 
year, being lower in fall and winter months [21]. Biophysical 
simulations using Niche Mapper also demonstrated that food 
intake and foraging rate were lower during these seasons. Lower 
T[, in winter compared to summer suggests that in winter lizards 
are under environmental constraints for thermoregulation and 
should operate as low-energy systems [20], [21]. Thus there 
appears to be a trade-off in choosing the use of terrestrial and/ or 
intertidal habitats by a lizard. In summer when T^. in the intertidal 
habitat is >25°C lizards use this microhabitat, extending their 
daily activity time and accepting the physiological cost of heat loss 
(i.e. a decrease in metabolic rate), while in seasons other than 
summer, the cost of heat loss appears to surpass the benefits and 



Table 4. ANOVA results for body size comparisons among sites, and sexes. 





Dependent variable 


Effect 


df 


MS 


F 


P-level 


SVL 


Site 


2 


0.00 


0.9 


0.34 




Sex 


1 


0.21 


100.9 


<0.001 




Site X Sex 


2 


0.00 


1.2 


0.31 




Error 


118 


0.00 








Site 


2 


0.14 


6.9 


<0.001 




Sex 


1 


2.57 


128.8 


<0.001 




Site X Sex 


2 


0.03 


1.7 


0.19 




Error 


126 


0.02 







SVL: Snout-vent length, Wt: body mass. 
doi:1 0.1 371/journal.pone.0097735.t004 
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lizards are restricted in their visits to the intertidal habitat. Thus 
the net ecological benefit establishes the use of difTerent 
microhabitats and consequently the activity period at seasonal 
(and geographical) scales. 

It is interesting to note that, although northern locations and dry 
conditions ofifer a considerable advantage in terms of time spent 
within optimal limits [44], the access of Kzards to wet habitats give 
them the advantage of more food, reducing their energy 
requirements. For these reasons, spatial and temporal differences 
in the access to the intertidal habitat could affect energy 
acquisition differentially for populations of Af. atacamensis. A 
latitudinal variation in the diet composition of A4. atacamensis was 
found, with an increase in the proportion of terrestrial items in a 
north to south gradient [18]. The diet of lizards from Medano 
(northernmost site) is composed almost exclusively of marine 
invertebrates and algae, thus lizards from this locality may be 
under a trade-off between feeding and thermoregulation, whii:h 
may in turn determine differences in the use of the thermal 
substrate. Due to the lower T^, except in the warmest months, 
lizards from Medano are restricted in visits to intertidal zones and 
this would reduce the ingestion of intertidal prey. This energy 
constraint might explain the lower body mass found in Medano 
lizards. Although longer seasonal activity is usually associated with 
larger body size [45], we found that hzards (both males and 
females) from Arrayan, and not from Medano, were larger. 

However, there are other possible contributing factors for the 
lack of relationship between activity duration and body size which 
should be considered. Together with a variation in temperature 
along the gradient there is a consistent increase in rainfall to the 
south, which promotes productivit}'. Lower abundance of tc'rres- 
trial invertebrates available for lizards may force northern lizards 
to exploit intertidal feeding because of the lack of insects in the 
extremely dry supralittoral zone. Thus even with more activity 
time they may be food limited, which could either stunt growth 
resulting in smaller body size or accelerate sexual maturation [46] , 
[47], which could also account for the observed smaller body sizes. 
Another hypothc-sis is that longer times at higher body temper- 
ature lead to higher metaboUc rates and therefore more food and 
water are needed to sustain those higher rates [48]. In fact, 
Maritza Sepulveda (unpublished data) found that lizards from 
Medano exhibited higher metabolic rates than lizards from 
Zenteno and Arrayan. If food and/ or ^vatcr are limited, greater 
metabolic expenditures in the face of limited food could also affect 
the potential growth. Finally, due to the lower temperature and a 
shorter activity season, lizards from Arrayan may need to extend 
the time of sexual development, resulting in larger body size at 
sexual maturity [13]. 

The activity pattern of M. atacamensis differs between juveniles 
and adults. We found that juveniles emerge before adults and 
become inactive earlier in the afternoon. At least three not 
mutually exclusive explanations can be hypothesized. First, 
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